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THE EFFECT OF A RING CURRENT ON THE BOUNDARY 
OF THE GEOMAGNETIC FIELD IN 
A STEADY SOLAR WIND 

By John R. Spreiter and Alberta Y. Alksne 


SUMMARY 


Approximate solutions are given for the shape of the boundary separating a 
steady neutral stream of ionized solar corpuscles from the combined magnetic 
fields of a three-dimensional dipole and an equatorial ring current. Results 
are presented for the traces of the boundary in the geomagnetic meridian plane 
containing the sun-earth line for several orientations of the latter relative to 
the dipole axis, and for the trace of the boundary in the geomagnetic equatorial 
plane for the case in which the dipole axis is normal to the sun-earth line. It 
is found that the presence of a ring current having values for the diameter and 
strength of the order proposed to explain the magnetometer data from Pioneer I 
and Pioneer V has the effect of greatly increasing the size, as well as altering 
the form, of the region within which the geomagnetic field is confined. 


INTRODUCTION 


The present paper reports the results of an extension of the theoretical 
study reported in references 1, 2, and 3 in which approximate results are deter- 
mined for the traces, in the geomagnetic equatorial plane and in the geomagnetic 
meridian plane containing the sun-earth line, of the cavity carved out of a 
steady neutral ionized solar corpuscular stream by interaction with a magnetic 
dipole representing the geomagnetic field. The novel feature of this extension 
is the inclusion of the effect of an equatorial ring current having properties 
similar to those of the model proposed by Smith, Coleman, Judge, and Sonett in 
reference k to represent the magnetometer data from Pioneer V and Explorer VI. 
These properties are that there exists, during quiet times, a westward flowing 
current of about 5xl0 6 amperes distributed over a large volume having the form 
of a toroidal ring 3 earth radii in cross-sectional radius with its center line 
situated in the geomagnetic equatorial plane at a distance of approximately 8 to 
10 earth radii. The magnetic moment of such a current system is of the same 
sign and order of magnitude as that of the main dipole field. Although the con- 
cept of a ring current is of long standing in the explanation of the decrease of 
the horizontal component of the magnetic field in the main phase of a magnetic 



storm, and values for the strength and radius similar to those stated above have 
also been deduced recently from cosmic-ray data by Kellogg and Winckler (ref. 5), 
magnetometer and other data from more recent space experiments with Explorer X 
(ref. 6 ) and Explorer XII (ref. 7 ) have failed to detect the presence of such 
substantial ring-current effects. As a result, the entire subject of the prop- 
erties and even the existence of a significant ring current must be regarded as 
an open question at the present time. 

It is evident that the presence of a ring current having a magnetic moment 
comparable to that of the permanent magnetic field of the earth should have the 
effect of greatly increasing the size, as well as altering the form, of the cav- 
ity. It is the purpose of this paper to present the results of a number of 
calculations undertaken to determine in a more quantitative manner the nature of 
the effects of such a ring current. A preliminary account of the present inves- 
tigation including plots of the results for the case in which the dipole axis is 
normal to the direction of the incident stream has been given in reference 8 . 
Those results are presented in greater detail in the present paper. Also pre- 
sented are additional results for the traces of the boundary in the above 
specified meridian plane for other orientations of the dipole axis relative to 
the stream direction. 

It should be recognized that the concept of a toroidal ring current with 
protons and electrons circulating round the geomagnetic axis at different speeds, 
although simple and of long standing, is, at best, not particularly precise. 
Singer has proposed a different concept in reference 9 in which the observed mag- 
netic variations and the associated ring current are explained in terms of the 
spiraling and drifting motion of trapped particles similar to those discovered 
shortly thereafter by satellite experiments by Van Allen and his colleagues (see, 
e.g., ref. 10). As is now familiar, the trapped particles spiral rapidly to and 
fro along the local magnetic lines between mirror points in the northern and 
southern hemispheres. At the same time the trapped particles drift round the 
earth, the protons to the west, the electrons to the east, thus setting up a 
westward electric current. This current is identified with the ring current. 
There is, in addition, a diamagnetic effect produced by the gyration of the par- 
ticles about the line of force. Superposition of these two effects gives the 
total magnetic influence, although Singer states in reference 11 that under cer- 
tain conditions, the diamagnetic effects are not important. 

It is assumed in the present investigation that the simple ring current 
model proposed by Smith, Coleman, Judge, and Sonett in reference k can be used to 
obtain an adequate representation of the basic geomagnetic field. It should be 
recognized that the values indicated for the strength and position of the ring 
current may be at considerable variance with those of the equivalent diamagnetic 
ring current which could also be used to represent the measured field, but which 
would be more consistent with the concept of the magnetic field arising from the 
motion of trapped particles. It is not necessary, however, for the determination 
of the terminal shape of the geomagnetic field to duplicate the properties of the 
ring current in all details. It is important, though, that a consistent system 
be employed so that the approximations introduced in proceeding from the measured 
values to an equivalent current system are removed when the process is repeated 
in reverse order to calculate the magnetic field. This has been done in the 
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present analysis, and all results presented herein are based on the consistent 
use of the model and associated numerical values proposed by Smith, Coleman, 
Judge, and Sonett. The only deviation is that the ring current of finite cross 
section is replaced by an idealized one of infinitesimal cross section. This 
change simplifies the computations considerably and should lead to little differ- 
ence in the results, provided the boundary of the cavity is farther than about 2 
or 3 earth radii from the idealized ring current. 


PRINCIPAL SYMBOLS 


a 

a e 

B 

B f 

Bp 

B s 

% 

Bp 

E 

i 

K 

k 

Mp 

Mi 

m 

n 

a 

n 

r, cp, e 

A 

r 


radius of ring current, cm (see fig. 1) 
radius of earth, cm 

total magnetic field, gauss 

magnetic field due to currents in boundary, gauss 
intensity of geomagnetic field at equator x 0.312 gauss 
intensity of field at the boundary, gauss 
B e r o 3 / M p ( see e<l. (9)) 

B r r o 3 / M p ( see e 1- (9)) 

complete elliptic integral of the second kind 
current, e.rn.u. 

complete elliptic integral of the first kind 

modulus of elliptic integrals 

dipole moment of earth = a 3 EL. 

e *0 

dipole moment of ring current = ira 2 i 
mass of proton - 1 . 67 xlCT 24 gin 
number of protons per cm 3 

unit vector in the direction of the outward normal to the boundary 
spherical coordinates (see fig. 1) 
unit vector in radial direction 
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unit of length defined by equation (8), 


cm 


v 

y, z 

xS yS z T 

a 

A 

9 

A 

M- 

P 


F 

i 

L 

N 

P 

R 

U 

9 

P 


velocity of plasma, cm/sec 

rectangular coordinates (geomagnetic, see fig, 1) 
rectangular coordinates (geographic) 
a/r 0 (see eq. (9)) 

unit vector in direction of increasing 6 

angle between direction of undisturbed plasma stream and the geo- 
magnetic equatorial plane (see fig, 1) 

Mq/Mp 


r/r 0 

angle between direction of undisturbed plasma stream and the out- 
ward normal to the boundary (see fig. l) 


Subscripts 


front 

pertaining to ring current 
lower 

neutral point 

pertaining to permanent, or dipole, field 

rear 

upper 

component in the direction of increasing 6 
component in the direction of increasing p 

FUNDAMENTAL ASSUMPTIONS AND EQUATIONS 


The basic concepts of the present study are classical and stem from a long 
series of investigations by Chapman, Ferraro, Dungey, and others (see refs. 12 and 
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13 for a r^sumd) undertaken to explain the connection between solar flares and 
geomagnetic storms. The fundamental assumption is that there exists a rarefied 
neutral ionized corpuscular stream, the solar wind, consisting principally of 
protons and electrons in equal numbers, which flows past the earth at hypersonic 
speeds. The particles are presumed to be of solar origin, and to be of uniform 
velocity in the undisturbed incident stream. The direction of the incident 
stream will be referred to for convenience of discussion as though it were coming 
directly from the sun. It should be understood, however, that the direction of 
the sun is actually immaterial in the analysis, and that only the direction of 
the undisturbed incident stream relative to the dipole axis is significant. 
Interaction between the solar wind and the permanent geomagnetic field is such 
that a cavity, bounded by a thin current sheath, is carved out of the stream. 

The solar wind is thus confined to the exterior, and the geomagnetic field to the 
interior, of the cavity. The latter region, dominated by the geomagnetic field, 
is now generally referred to as the magnetosphere. 

The particles of the solar wind are considered to move in straight lines up 
to the boundary of the cavity where they are, in effect, specularly reflected and 
returned to the stream with a direction of motion different from that which they 
possessed in the incident stream. By this process they exert a pressure 
2mnv 2 cos 2 ^ on any element of the boundary for which cos < 0. The quanti- 
ties m, n, and v represent the mass, number density, and velocity of the pro- 
tons in the undisturbed incident stream, and ^ represents the angle between the 
free-stream velocity vector and an outward normal to the surface. An element of 
the surface that fails to coirply with the condition that cos < 0 is shielded 
from the stream and experiences no pressure. 

Dungey (ref. 13 ) has investigated the conditions that prevail in the current 
sheath bounding the cavity and has shown that the particle pressure is balanced 
by the magnetic pressure B s 2 /8jt where B s is the total (tangential) magnetic 
field at the surface of the cavity. These considerations lead to the following 
relation which must be satisfied at the boundary: 

B s 2 /8tx = 2mnv 2 cos 2 t (l) 

with m, n, and v expressed in c.g.s. units and B s expressed in gauss. 

The total magnetic field B in the cavity is considered to be the sum of 
the permanent magnetic field Bp of terrestrial origin, the field Bj_ induced 
by the ring current, and the field B' due to the currents in the boundary. It 
is a fundamental assumption of the theory that the total field inside the cavity 
satisfy the magnetic field equations for a vacuum 

div B = 6 , curl B = 0 ( 2 ) 

together with the boundary conditions that the normal component of B vanish and 
the tangential component of B be given by equation (l) at the surface of the 
cavity. It is also necessary that the solution possess the appropriate singular- 
ities in the cavity that are required to represent Bp and Bj_ . The first of 
these will be considered to be represented by a single magnetic dipole, thus 
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(3) 


Bp = + rBp = -(M p /r 3 )(© sin 0 + r2 cos 0) 

where the coordinate system is fixed with respect to the dipole as illustrated in 
figure 1, 0 and r are unit vectors in the direction of increasing 0 and r, and 
the magnetic moment of the dipole is given by Mp = a e 3 Bp where ae represents 
the radius of the earth (6.37x10 s cm) and B-,-, represents the intensity of the 

geomagnetic field at the magnetic equator (0.312 gauss). The second is given by 
(see, for instance, ref. 1^) 



+ rBi 


r 


2Mq 


:rta 2 (a 2 + r 2 + 2ar sin 



t fjr a g + r g 
sin 0 \ a 2 + r 2 


2a g sin g 0 
2ar sin 0 



A 

- r 



2a 2 cos Q 
+ r 2 - 2ar sin 6 



where E and K are complete elliptic integrals having the modulus 

1/2 

f Aar sin 9 

k = 


+ r 2 + 2ar sin Q. 




(5) 


and a and Mj_ represent the radius and magnetic moment of the ring current. The 
latter is defined by 

Mi = ira 2 i (°) 

where i represents the current measured in e.m.u. and considered positive when 
flowing westward around the earth. One e.m.u. of current is equal to 10 amperes. 

It is a property of the above equations that B cannot vanish over any 
region of finite extent in the interior of the cavity. It follows that the 
boundary of the cavity must be of such form that cos i < 0 everywhere, with 
zero values occurring only at isolated points. This condition together with the 
above equations suffices to specify completely the form for the boundary of the 
cavity and the properties of the magnetic field contained therein. 


DERIVATION OF AN APPROXIMATE DIFFERENTIAL EQUATION FOR THE 
COORDINATES OF THE BOUNDARY OF THE CAVITY 


Although the simplified physical model described above leads to a completely 
defined mathematical problem, the exact solution has still to be determined, even 
for the case of the dipole alone. The present analysis, therefore, is based on 
the use of the approximation introduced and commented upon recently by Beard 
(ref. 15), Ferraro (ref. l6), and Davis and Beard (ref. 17). It is judged on the 
basis of previous experience with the case of the dipole alone that results 
obtained in this way display nearly all the essential features of the exact 
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solution, although the coordinates of the boundary may be in error by a few- 
percent one way or the other at various points. This statement is based princi- 
pally on comparisons of exact and approximate solutions for the analogous problem 
in two dimensions given by Hurley (ref. l 8 ) , and by Spreiter and Briggs (refs. 1 
and 2 ) . 

The essential concept that leads to the simplif ication achieved by Beard is 
that of relinquishing the condition that the normal component of B vanish at 
the boundary of the cavity and replacing it by the approximate relation that B s 
be equal to twice the tangential component of the permanent magnetic field B^ 
at the same point. Ferraro subsequently suggested that it would be better to 
replace the factor 2 by 2 f where f is a constant, the value for which is to be 
determined at the end of the calculation by matching conditions at some partic- 
ularly significant point. These concepts have been carried over to the present 
study, for which they lead to the approximate relation 

B g 2 = ^f 2 [n X (Bp + Bp ] 2 (7) 

where n is a unit vector in the direction of the outward normal to the surface 
of the cavity. Ferraro proceeds to present a simple illustrative two-dimensional 
example involving flow past a current bearing wire, and shows that a reasonable 
procedure for the estimation of f leads to the value 0.68 for that particular 
case. It was subsequently shown elsewhere that the value for f increases to 
about 0.91 when flow past a two-dimensional dipole is considered (refs. 1 , 2 , and 
l 8 ), and to 1.29 for a three-dimensional dipole immersed in a stationary plasma 
(ref. 19)* Since the linear dimensions of the cavity are proportional to the 
first power of f for the wire problem, to the square root of f for the two- 
dimensional problem, and to the cube root of f for the three-dimensional dipole 
problem, the use of these values instead of unity leads to reductions of all 
dimensions of the cavities around the wire and the two-dimensional dipole by 
32 percent and 5 percent, respectively, and to an increase of the dimensions of 
the cavity around the three-dimensional dipole by 9 percent. Corresponding 
results have not been determined for the present case, however, and f is 
equated to unity in all numerical examples presented herein. 

It should be observed that the expression for B s given by equation ( 7 ) is 
as indicated by Davis and Beard in reference 17 and, due to the correction of an 
error, differs from that given previously in references 1, 2, 8, 15 , and else- 
where. As noted by Davis and Beard, the difference disappears when attention is 
confined to the meridian plane containing the sun-earth line, and to the equa- 
torial plane for the case in which the dipole axis is normal to the direction of 
the incident stream. Solutions given previously for these planes are thus cor- 
rect in spite of their derivation, but results for the coordinates of the remain- 
der of the surface must be recalculated in order to be consistent with 
equation (7)- 

Equations ( 3 ) through (7) suffice to provide a relation expressing the left- 
hand side of equation (l) in terms of the coordinates of the surface of the cav- 
ity, the dipole moment of the permanent magnetic field, and the radius and mag- 
netic moment of the ring current. The variable part, cos of the right side of 
equation (l) can also be expressed in terms of the coordinates of the surface. 



Combining these relations leads to a partial differential equation for the radial 
coordinate r of the surface expressed as a function of the angular coordinates 
Q and cp. The equation is lengthy, however, and it is desirable before display- 
ing it to achieve a certain economy resulting from a reduction in the number of 
parameters by introducing the following quantity for the unit of length 


r 


o 



a c 


4f%o 2 

•l6rtmnv 2 . 


1/6 


( 8 ) 


together with the following ratios 


P = r/r 0 , a = a/r c , n = Mj/Mp, B 0 = B 0 r o 3 /Mp, B p = B r r 0 3 /M r 


(9) 


The governing equation can thus be written 


* Bi p) + (% " Bl e) /a 0 ' 2 


p 2 sin 2 9 




+ B 

e. 


v\VA/3£ 

' p a© 


sin cp sin 0 cos A - cos 0 sin A - ■ (sin cp cos 0 cos A + sin 9 sin A) 

p 3e 


dp 


p sin 9 dcp 


cos cp cos A 


= cos 2 t 


where 


B. 


^9 


B-n = -sin 
P 0 


2|i 


e/p 3 , b p = -2 cos e/p : 


( 10 ) 


(ii) 


ira 2 (a 2 + p 2 + 2ap sin 9)^'^ sin 9 

-4pa^E cos 9 


i/a 


B-i 


K a 2 + P 2 - 2a 2 sin 2 9 g 
a 2 + p 2 - 2ap sin 9 


> (12) 


P 2/2 2 


5ta 


(a^ + p^ + 2ap sin e) 1/2 ( a 2 + p 2 - 2ap sin 9) 


where E and K are complete elliptic integrals having the modulus 

N l/2 

k _ ( 4ap sin 9 


a 2 + p 2 + 2ap sin 9. 


(13) 
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where A denotes the angle in the yz plane between the direction of the inci- 
dent stream and the normal to the axis of the dipole as shown in figure 1. The 
quantity r Q bas physical significance as the geocentric distance along the sun- 
earth line to the boundary of the geomagnetic field for the case i = 0, A = 0. 

The dimensionless radius vector p defining the form of the boundary of the cav- 
ity as a function of 0 and cp thus depends on the parameters and A. The 
determination of the size of the cavity for any particular combination of density 
and velocity of the incident stream requires, in addition, the calculation of the 
unit of length r Q - 

When Sp/ckp = 0, as it is in the meridian plane containing the sun-earth 
line, the first term on the left side of equation (10) is zero and it is possible 
to take the square root of both sides of the equation. In so doing, it is impor- 
tant that the ± sign is used in such a way that both sides of the equation are 
of the same sign. The term within the square bracket on the right side is equal 
to -cos ^ and is therefore positive. The term within the square bracket on the 
left side may be of either sign, however, depending upon whether the tangential 
component of Bp + Bq is oriented in the direction of increasing or decreasing 
0. Both cases occur in the desired solution, since all field lines on the sur- 
face of the cavity diverge from and converge to a pair of isolated neutral points. 
These points are of additional interest because the magnetic field lines turn 
abruptly there and extend to the earth. The intersections of these lines and the 
earth T s surface define a pair of isolated points in the vicinity of the geomag- 
netic poles through which pass all of the field lines that lie in and define the 
boundary of the cavity. These points are of special significance because they 
locate the geographical areas into which can precipitate charged particles 
initially trapped in the vicinity of the boundary of the geomagnetic field. 


SOLUTION FOR THE MERIDIAN PLANE CONTAINING DIPOLE AXIS AND SUN-EARTH LINE 


The solution of equations (10) through ( 13 ) is unquestionably complex, but 
can be simplified considerably if attention is confined to determining the trace 
of the boundary of the cavity in the meridian plane containing the dipole axis 
and the sun-earth line, that is, the plane along which cp = ±jt/ 2. Along this 
plane Sp/Scp vanishes by reason of symmetry, and the governing equation reduces 
to an ordinary differential equation in which bp/bd can be replaced by dp/d0. 
This equation can be solved for dp/d0 to yield 


dp = 


P 


d0 


(sin 0) [cos A sin (±;rt/2) + 
(cos 0) [ cos A sin (±:rt/2) ± 


(l/p 3 )] - cos 0 sin A ± Bj_ 

0 

(2/p 3 )] + sin Q sin A + B ip 


( 14 ) 


where the sense of the plus and minus signs has been retained to be the same as 
in the preceding equations by writing sin (±tt/2) for sin (cp = ±jt/2) = ±1 in 
the right-hand member. It is convenient to consider the two rearranged families 
of solutions associated with the following equation, equivalent to equation (l4) 
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sin [0 - A sin (±jt/ 2 )] + (sin 0 )/p 3 ± Bj~ 

-JL (15) 

cos [0 - A sin (±rt/ 2 )] ± 2 (cos 0 )/p 3 + B^p 

^ 4 

If the intensity of the ring current is zero, Bi^ = Bj_p = 0, and the solution 

for any A can be determined in analytic form as shown in references 1 and 2 . It 
was found that the upper signs give a result appropriate for the nose or front 
portion of the boundary, and the lower signs a result appropriate for the rear 
portion. The results become particularly simple for the case in which the inci- 
dent stream is normal to the dipole axis so that A = 0. In this case the coor- 
dinates of the front portion of the boundary are given by 

p = l ( 16 ) 

and those of the rear portion by 




The curves representing the front and rear portions of the boundary join together 
at an angle given for the upper half plane by 



— radian * 19 - 1 ° 

3 


( 18 ) 


This point and the corresponding point in the lower half plane are of particular 
physical significance because they correspond to the neutral points described in 
the preceding section. 

The greater complexity inherent in equation (15) when the intensity of the 
ring current is not zero precludes solution in analytic form, but results can 
still be obtained by standard numerical techniques. There remains, however, the 
discussion of the choice of the appropriate integration constant or, what is 
equivalent, the appropriate combination of p and 0 to use as starting values in 
the numerical integration. The appropriate choice for the starting values 
required to define uniquely the desired solution that satisfies all the auxiliary 
conditions of the problem is not immediately evident, since there is no point on 
the boundary for which the coordinates p and 0 are known from a priori consid- 
erations. Examination of the properties of the integral curves for equation (15) 
discloses, however, that only one of the many alternative solutions satisfies the 
condition that the cavity extends a finite distance from the ring current in the 
direction of the sun, and that the lateral dimensions increase steadily with dis- 
tance downstream from the apex so that cos < 0 at all points. 

This statement is illustrated most readily by consideration of a case in 
which A = 0. Plots of the integral curves for the upper half plane are shown in 
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figure 2 for the case in which a = 1 and p = 1. The corresponding curves for 
the lower half plane need not be considered for the case of A = 0, since the 
form of the hollow is symmetric about the sun-earth line. Although the integral 
curves shown in figure 2 are for a particular pair of values for a and p, they 
are typical because integral curves for other values of these parameters are 
qualitatively similar. Also indicated on these plots are the positions of the 
three singular points of equation ( 15 ) that are in or adjacent to the upper half 
plane. These points are designated according to their general location by the 
terms, front Sf, rear Sr, and upper Sjj- Their coordinates are determined by 
the condition that the numerator and denominator of the right-hand members of 
equation (15) vanish simultaneously. Their significance derives from the fact 
that the only integral curves that satisfy all the stated requirements for an 
acceptable solution pass through a singular point. Thus, the only integral curve 
that intersects the sun-earth line at a finite and nonvanishing distance upstream 
from the ring current is that shown on part (a) of figure 2, which passes through 
Sf- This curve is, therefore, the only one that can represent the trace of the 
cavity boundary in the vicinity of the apex. This curve cannot represent the 
trace of the boundary for cp = -jt/2, however, since it turns away from the direc- 
tion of the corpuscular stream over the pole and fails to satisfy the condition 
that cos \|r < 0. The integral curves shown on figure 2(b) reveal that there is 
one and only one integral curve that can be joined to the curve which represents 
the front portion of the boundary, and that extends to infinity in the rearward 
direction with cos \|r < 0 at all points. It is the integral curve shown in 
part (b) that passes through Su- 

A plot of the resulting form for the trace of the cavity boundary is shown 
in figure 3 for the case in which a = 1 and p = 1. Also indicated on this plot, 
by dashed lines, are the corresponding results defined by equations ( l6 ) and (17) 
for the case in which there is no ring current. 

If A / 0, the symmetry about the equatorial plane disappears and it is 
necessary to determine integral curves for the lower half plane as well as the 
upper. A representative set of integral curves for A = 11-5° is shown in fig- 
ure 4 for the case a = 1, p = 1. Solutions for small A are determined in a 
manner similar to that just described for A = 0 by selecting the exterior por- 
tions of the integral curves that pass through the singular points. It is thus 
necessary to consider the lower singular point Sp which is situated at the same 
geocentric distance as Sp but in the opposite direction from the origin as 

indicated in figure 4. The resulting curve for the trace of the boundary of the 

geomagnetic field in the meridian plane containing the sun-earth line is shown in 
figure 5 for the case a = 1, p = 1. This curve represents the only solution of 
equation (15) for these values for a and p that is continuous and satisfies the 
condition cos < 0 at all points. As in figure 3* the dashed line indicates 
the results from references 1, 2, and 3 for the case in which there is no ring 

current. It is of interest to observe that the angular coordinates of the neu- 

tral points, at which the upper and lower curves join the curve from the forward 
portion of the boundary, change at a slower rate than the angle between the 
stream direction and the dipole axis. 

As A is increased it reaches a critical value A cr , dependent upon the 
values for a and p, at which the upper singular point coincides with the upper 
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neutral point. If A is increased beyond A cr the procedures described above 
to determine the upper rear portion of the boundary no longer apply because the 
condition cos ^ < 0 is violated on the portion of the boundary immediately 
upstream of the indicated neutral point. The analogous situation for the case of 
the dipole alone was encountered in reference 1 when A exceeded 35-6°. A typi- 
cal example in which a ring current is present is provided by considering the 
case a = 1, p = 1, A = 34. 5° • The integral curves for this case are shown in 
figure 6 . The dashed line on figure 6 marks the locus of the points at which the 
integral curves are parallel to the stream, that is, where cos i|r = 0. For the 
case of the dipole alone this would occur along a straight line, Q = Qq^j for 
cp = Rjt/ 2, and 0 = 0s^ for cp = -jt/2. The addition of the ring current alters 
the situation and the locus is now given by 

sin (e - \) gf— - B ip ^) + cos (6 - 7 \) / - B l0 ) = 0 ( 19 ) 

This equation applies to both families of curves and to both the cp = +rc/2 and 
the cp = -tc/ 2 planes. The only integral curves for which cos ^ / 0 at the 
intersection with this line are the four members of the family shown in fig- 
ure 6 (b) which pass through the singular points Su and Sp. As in the case of 
the dipole alone, the desired solution can be obtained by joining the appropriate 
curves representing the forward and rearward parts of the boundary at the point 
of intersection with the line given by equation (19) • Both portions of the 
boundary trace are parallel to the direction of the incident stream at this point 
and satisfy the condition that cos ^ < 0 everywhere. 

The resulting solution for the trace of the boundary in the meridian plane 
is shown in figure 7 . As in figures 3 and 5, the singular and neutral points are 
represented by circles and squares and the dashed lines show the corresponding 
results for the case in which the ring current is absent. 

Figures 3, 5, and 7 show the effect of a ring current on the size and shape 
of the boundary for three different orientations of the dipole axis relative to 
the stream direction. The principal effect is an over-all enlargement, with some 
distortion near the front, where the boundary is closest to the ring current. 


SOLUTION FOR EQUATORIAL PLANE FOR A = 0 


For the special case in which A = 0, simplification comparable to that in 
the meridian plane occurs if attention is confined to determining the trace of 
the boundary in the equatorial plane. Along this plane 0 = tt/ 2, 6p/d0 vanishes 
by reason of symmetry, and the governing equation reduces to an ordinary differ- 
ential equation in which dp/6 cp can be replaced by dp/dcp. This equation can be 
solved for dp/dcp, and the result is as follows if attention is restricted to the 
interval tt/2 < cp < 3 tt/2 


12 





dp 

drp 


p 6 sin cp cos cp + (l - b iP 3 )Vp s - (d - b iP 3 
~ 6 cos 2 cp - (l - bj_p 3 ) 2 


2 -» 


( 20 ) 


where 


, _ 2d ( K 

1 *a a W P 



= (BiJ 


0= It/ 2 


( 21 ) 


in which E and K are complete elliptic integrals having the modulus 

n , , 1/2 

k = 2(a P ) (22) 

a + p 

If the intensity of the ring current is zero,, q and, hence, bj_ vanish, and equa- 
tion (20) reduces to the form given in references 1 , 2, and 3 for the dipole 
alone . 

The integration constant or starting values to be used in the integration of 
equation (20) must be such that the curve describing the desired solution inter- 
sects the sun-earth line at the same point as the curve describing the trace of 
the boundary of the hollow in the meridian plane. This signifies that the inte- 
gral curve of equation (20) that is to be selected to represent the boundary of 
the hollow must pass through the front singular point of equation ( 15 ) with upper 
signs. That this is the only integral curve of equation (20) which intersects 
the sun-earth line at a finite distance from the dipole singularity, and hence is 
an acceptable solution of the original physical problem, is apparent from consid- 
eration of the properties of a set of integral curves of equation (20) illus- 
trated in figure 8 for the case defined by a = 1 and \± = 1. 

Selection of the particular integral curve described in the preceding para- 
graph yields the curve defining the trace of the cavity boundary in the equato- 
rial plane for the case in which A = 0. An isometric view of this curve is 
illustrated in figure 9 > together with the corresponding curve from figure 3 for 
the trace of the boundary in the meridian plane containing the sun-earth line, 
for the special case a = 1, p = 1. The corresponding curves for the case of the 
dipole alone are also indicated in figure 9 by the dashed lines. 


RESULTS AND DISCUSSION 


The procedures described in the preceding sections enable the calculation of 
traces of the boundary of the geomagnetic field in the equatorial plane for the 
case A = 0, and in the meridian plane containing the dipole axis and the sun- 
earth line for various A, for any given set of values for the velocity and num- 
ber density of the solar wind and for the strength and radius of the ring 
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current. The selection of appropriate values to achieve realistic results 
involves some uncertainties , however, not only because knowledge is meager, but 
also because the magnitudes of the quantities are believed to vary considerably 
with time. The results that follow are presented, therefore, for cases involv- 
ing what are believed to be representative values for these quantities together 
with results for other cases having smaller and larger values of the governing 
parameters . 

Smith, Coleman, Judge, and Sonett have shown in reference 4 that the magne- 
tometer data from Explorer VI display variations that can be matched closely by 
results computed assuming the presence of an equatorial ring current of 5 million 
amperes located at a distance of 60,100 km from the center of the earth and dis- 
tributed over a cross section having a radius of 3 earth radii or less. This 
finding forms the basis for the computed forms for the boundary of the cavity 
that are presented in figures 10 through 13. The results shown in these figures 
and tabulated in tables I and Ilia are for a ring current of infinitesimal cross 
section located 60,000 km from the center of the earth. 

In the results presented in figure 10 , the solar corpuscular stream is 
assumed to have a number density n in protons/cm 3 and velocity v in cm/sec 
so that the product nv 2 is equal to 10 1S . This value would be possessed, for 
instance, by a stream having a number density of 4 protons/cm 3 and a velocity 
of 500 km/sec. Such values are representative of conditions observed by 
^Mariner II during quiet times (ref. 20 ). Results are shown for three different 
values for the current, namely, 2 , 5 ^ and 10 million amperes. The results for no 
ring current are also indicated by a dashed line. Angles of incidence extending 
from 0 ° to 34.5° are included so as to cover the range of angles which the geo- 
magnetic equatorial plane makes with the sun-earth line during the course of the 
year. The results for A = 0 are also shown in isometric projection in fig- 
ure 11 in order to help visualize the three-dimensional shape of the boundary. 

The corresponding plots for other A are not shown because results are available 
for only the meridian plane. The substantial effects of the ring current in 
enlarging the size of the cavity are clearly illustrated by the results. It can 
be seen, for example, that when A = 0 the distance to the termination of the 
geomagnetic field along the sun-earth line increases from about 8.8 to 12.4 earth 
radii with the addition of a ring current of 5 million amperes to the basic 
dipole field. 

The number density and velocity of the solar wind also have a considerable 
effect on the dimensions of the cavity. The magnitude of the effect is illus- 
trated in figures 12 and 13 for the case in which the strength of the ring cur- 
rent is 5 million amperes and the product nv 2 takes on several different values 
between 10 15 and 10 17 with n in protons/cm 3 and v in cm/sec. As anticipated, 
the size of the cavity diminishes rapidly as the intensity of the solar wind 
increases. Results are not given for nv 2 greater than 10 17 because part of the 
boundary approaches so near to the idealized ring current that the simplifying 
assumption of infinitesimal radius becomes totally inadequate. 

Smith, Coleman, Judge, and Sonett have also analyzed the magnetometer data 
from Pioneer V, which was launched during the recovery phase of a magnetic storm, 
and have shown in reference 4 that it could be duplicated by results computed 



assuming the presence of an equatorial ring current of 5 million amperes located 
at a distance of 50,000 km from the center of the earth and distributed over a 
circular cross section having a radius not less than 3 earth radii. These values, 
which differ from those indicated by the data from Explorer VI only in the dis- 
tance from the earth to the ring current, form the basis for the selection of the 
parameters defining the results for the cavity boundary shown in figures Ik 
through 17 and tabulated in tables II and Illb. The results are for a ring cur- 
rent of infinitesimal cross section located 50*000 km from the center of the 
earth. The other parameters are the same as those for which results are shown in 
figures 10 through 13. Thus, the results presented in figures l 4 and 15 illus- 
trate the effect of varying the strength of the ring current, with the number 
density and velocity of the solar corpuscular stream held fixed at such values 
that nv 2 is equal to 10 16 . Those presented in figures l6 and 17 illustrate the 
effect of varying the intensity of the corpuscular stream with the strength of 
the ring current maintained fixed at 5 million amperes. It is evident, as is to 
be expected, that the dimensions of the cavity diminish somewhat as the distance 
from the earth to the ring current is reduced, other quantities remaining the 
same. The effect is not overwhelming, however, and the results for the two cases 
are similar in most qualitative features. 

Since the coordinates of the neutral points are of considerable interest, 
and they cannot be determined accurately from table I or II, a separate listing 
has been prepared and is presented herein as table IV. 

All equations and results presented in the preceding part of this paper are 
given from the point of view of an observer fixed with respect to the geomagnetic 
coordinates. It is, however, of equal interest to examine the results from the 
point of view of an observer fixed with respect to the direction of the solar 
wind. Figure l8 shows results assembled in this way for the case in which the 
stream direction is normal to the geographic axis and the geomagnetic axis makes 
an angle therefrom of ±11.5°, as occurs during the course of the day at the equi- 
noxes, provided the stream direction is assumed parallel to the sun-earth line. 
Corresponding results for the case in which the geographic equatorial plane is 
inclined 23° from the stream direction, such as occurs at the solstices, are 
shown in figure 19. These results show that the boundary of the geomagnetic 
field is only slightly affected by the diurnal wobbling of the geomagnetic dipole 
axis relative to the direction of the solar wind. This is particularly true for 
the dipole alone. The results for this case show that the only part of the 
boundary that has a marked effect of orientation is that in the Immediate vicinity 
of the neutral points. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Mar. 28 , 1963 
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TABLE I . - COORDINATES OF THE TRACE OF THE BOUNDARY OF THE GEOMAGNETIC FIELD IN 
THE MERIDIAN PLANE CONTAINING THE DIPOLE AXIS AND THE SUN-EARTH LINE FOR A 
RING CURRENT LOCATED AT A GEOCENTRIC DISTANCE OF 60,000 KM 


| (a) Stream 

= 10 ia protons /cm 

sec 2 ; ring current 

= 2 million amperes 

| (b) Stream = 2 - 5x10 16 protons/cm 

sec 2 ; ring current = 2 million amperes | 

\ X 

0 

11.5 

23.0 

34.5 

0 

U--5 

23.O 

34.5 

0 

r* 5 

1 23-0 

3^-5 

0 

11.5 

23.0 

34.5 

deg\ 

<p = +rr/2 


9 = - 

*/2 


I <p = +ir/2 


cp = - 

rt/2 


0 

11.55 

12.87 

14.65 

17-26 

11-55 

12.87 

14.65 

17.26 

9-84 

10.96 

12.49 

l4.Bi 

9.81T 

IO.96 

12.49 

14.81 

5 

10.86 

12.01 

13.52 

15.72 

12.29 

13.83 

15.94 

19.09 

9.26 

10.23 

11-52 

13.50 

10.47 

II.78 

13.59 

16.37 

10 

10.23 

11.22 

12.52 

14-39 

13.12 

14.91 

17.44 

21.30 

8.72 

9.56 

10.67 

12.40 

11.17 

12.71 

14.86 

18.26 

15 

9.65 

10.51 

U.62 

13-26 

14-03 

16. 14 

19.19 

24.03 

8.22 

8.95 

9.90 

11.47 

11.95 

13.76 

16.36 

20-59 

20 

9.29 

9-98 

IO.89 

12-35 

15.06 

17.57 

21.30 

27.50 

7.92 

8.59 

9-46 

IO.76 

12.83 

14.97 

18.16 

23-55 

25 

9-31 

9-95 

10.79 

11.90 

16.23 

19.24 

23-88 

32.07 

7-94 

8.57 

9.38 

10.43 

13.83 

16.39 

20.36 

27.47 

30 

9.34 

9.94 

10.72 

11.73 

17.59 

21.23 

£7. 13 

38.39 

7-97 

8.56 

9-32 

10.28 

14.96 

18.10 

23-14 

32.88 

35 

9.38 

9.94 

10.67 

11.59 

19.16 

23.66 

31.36 

47.76 

8.00 

8.57 

9.28 

10 -16 

16.33 

20.17 

26.75 

40-91 

Uo 

9-43 

9.96 

10.63 

11.47 

21-04 

26.69 

37-14 


8.04 

8.59 

9-25 

10.07 

17.93 

22.76 

31.69 


45 

9-50 

9-99 

10.61 

11.38 

23.31 

30.61 

45-50 


8.10 

8.62 

9-25 

10.00 

19.87 

26.11 

38.83 


50 

9-58 

10-04 

10-62 

11.31 

26.13 

35.87 



8.17 

8.68 

9-27 

9.96 

22.28 

30.61 



55 

9-69 

10.12 

10.65 

11.26 

29-73 

43-35 



8.27 

8.76 

9-32 

9.95 

25.36 

37-00 



6o 

9.83 

10.23 

10.70 

11.27 

34.51 




8.41 

8.88 

9.40 

9.98 

29.44 




65 

10.02 

IO.38 

10-79 

11.28 

41.18 




8.61 

9.06 

9-53 

10-04 

35-15 




70 

10.26 

10.56 

10-92 

11.33 

51.19 




8.90 

9.30 

9-72 

10.15 

43.70 




75 

10.56 

10.79 

11.07 

11-40 





9-31 

9.64 

9.96 

IO.30 





8o 

10.87 

11.03 

11.22 

11.46 




41.20 

9-84 

10.04 

10.25 

10.48 




35-18 

85 

11.12 

11.20 

11.32 

11.47 



56.21 

34.64 

IO.31 

10.39 

10A9 

10.61 



47.98 

29-58 

90 

11.22 

11.23 

11.29 

11.38 



44.39 

30.00 

10.49 

10-50 

10.52 

10.57 



37.89 

25.61 

95 

11.12 

11.08 

11.07 

11.10 



36.82 

26.54 

IO.31 

10.24 

10.19 

10.15 



31.43 

22.65 

100 

IO.87 

IO.76 

10.68 

10.65 


47.71 

31-56 

23.85 

9.84 

9.63 

9. 41 

9-14 


40.73 

26-93 

20.35 

105 

10.56 

10-37 

10.22 

10.12 


39.00 

27.68 

21.69 

9-31 

8.97 

8-57 

8.07 


33-29 

23.62 

18.49 

110 

10.26 

10.01 

9.80 

9.67 

51-19 

33-07 

24.69 

19-90 

8.90 

8.49 

8.07 

7-64 

43.70 

28.22 

21.06 

16.95 

115 

10-02 

9.72 

9.49 

9-33 

41.18 

28-75 

22.31 

18.40 

8.61 

8.19 

7-80 

7-45 

35.15 

24.52 

19.02 

15-66 

120 

9.83 

9-51 

9.26 

9.10 

34.51 

25-46 

20.36 

17.11 

8.4l 

7.99 

7-63 

7-34 

29.44 

21.71 

17-34 

14-54 

125 

9-69 

9.34 

9.08 

8.91 

29.73 

22.86 

18-73 

15.98 

8.27 

7.86 

7.51 

7-25 

25.36 

19-48 

15-94 

13-58 

130 

9-58 

9-21 

8.94 

8.77 

26.13 

20.74 

17.34 

15-00 

8.17 

7-75 

7-42 

7.16 

22.28 

17-67 

14-75 

12.75 

135 

9.50 

9-11 

8.82 

8.64 

23.31 

18-99 

16. 14 

14.13 

8.10 

7.67 

7-3^ 

7.H 

19.87 

16.17 

13-73 

12-01 

110 

9-43 

9-02 

8.72 

8.53 

21.04 

17-50 

15.09 

13.36 

8.04 

7.61 

7.28 

7-05 

17.93 

14-90 

12.84 

11-35 

145 

9.38 

8.95 

8.63 

8.43 

19-16 

16.23 

14.17 

12-66 

8.00 

7-55 

7-22 

6.99 

16.33 

13.82 

12.06 

10-77 

150 

9.34 

8.88 

8-55 

6-33 

17-59 

15-12 

13-35 

12.03 

7-97 

7-50 

7.16 

6.92 

14.98 

12. &7 

II.36 

10.24 

155 

9-31 

8.83 

8.47 

8.24 

16.23 

14.14 

12-61 

11.46 

7-94 

7-46 

7-10 

6.86 

13.83 

12-04 

10.74 

9-75 

160 

9-29 

8.77 

8.U) 

8.14 

15.06 

13-27 

11.94 

10-93 

7-92 

7-42 

7-05 

6-79 

12.83 

11.31 

10.17 

9-31 

165 

9-65 

8-97 

8.44 

8.05 

14.03 

12.49 

11-33 

10-44 

8.22 

7-64 

7-19 

6.82 

11-95 

10.64 

9.65 

8.89 

170 

IO.23 

9-46 

8.85 

8.37 

13-12 

11-79 

10.77 

9-90 

8.72 

8.06 

7-55 

7.13 

11-17 

10.04 

9-18 

8.50 

175 

10. e6 

9-98 

9.29 

8-75 

12.29 

11.14 

10-25 

9-55 

9-26 

8.51 

7-92 

7-45 

10.47 

9-49 

8-73 

8.14 

180 

11.55 

10-54 

9.76 

9-14 

U-55 

10-54 

9.76 

9-14 

9-84 

8.98 

8.31 

7-79 

9-84 

8.98 

8.31 

7-79 


(c) Stream = 10 1S 

protons/cm sec 2 ; ring current = 5 million amperes 

(d) Stream = 10 le 

protons/cm sec 2 ; 

ring current = 5 million amperes J 

0 

18.62 

20.79 

23.73 

27.90 

18.62 

20.79 

23-73 

27-90 

12.25 

13-73 

15.73 

18.74 

12.25 

13-73 

15-73 

18.74 

5 

17.52 

19-41 

21.91 

25-40 

19.82 

22.35 

25-81 

30. 87 

11.53 

12.81 

14.52 

17-08 

13-04 

14.76 

17-12 

20.71 

10 

16.51 

18.15 

20-30 

23-25 

21-16 

24.10 

20.23 

34.45 

IO.87 

11-99 

13-46 

15-68 

13-93 

15-92 

18.73 

23-11 

15 

15-59 

17-02 

18.86 

21.39 

22.65 

26.10 

31-07 

38.87 

10.26 

11.24 

12-51 

14.50 

14-91 

17.25 

20.63 

26.07 

20 

15-03 

16-01 

17-57 

19.79 

24-33 

28.41 

34.48 

44.47 

9.99 

10.79 

11.84 

13.56 

16.03 

18.80 

22.92 

29.85 

25 

15-08 

15-99 

17-21 

18.86 

26.24 

31.12 

36.66 

51.86 

10.04 

10.79 

II.76 

13.03 

17.31 

20.61 

25-72 

34.82 

30 

15.16 

15-99 

17.10 

18.59 

28.44 

34-35 

43.92 

62.09 

10.11 

10.81 

11.70 

12-86 

I8.78 

22.78 

29.25 

4i.yi 

35 

15-24 

16.00 

17-02 

10-37 

31.00 

38.28 

50.78 

77-23 

10.20 

10.85 

II.67 

12.73 

20.51 

25-42 

33-85 

51.90 

10 

15-35 

16. 04 

16.96 

18.18 

34.04 

43-19 

60.12 

; 

10-31 

IO.91 

IV 67 

12.64 

22.56 

23-72 

*40.11 


45 

15-46 

16.09 

16.92 

18.02 

37-71 

49-52 

73-65 


10-44 

11.00 

0.70 

12.57 

25-03 

32-97 

49.17 


50 

15.60 

16.15 

16.90 

17.89 

42.27 

53-03 



10-60 

11.12 

11.76 

12.54 

28.10 

38.67 



55 

15.74 

16.23 

16.89 

17.76 

U8.09 

70.12 



10.80 

11-27 

11.84 

12.54 

32.01 

46-77 



60 

15.89 

16.31 

16.90 

17-69 

55-31 

88-74 



11.04 

11.45 

II.96 

12.57 

37-19 

59-21 



65 

l6.04 

16. 40 

16.91 

17.61 

66.60 




11-31 

11.66 

12-09 

12.63 

44.41 




70 

16.19 

16.48 

16.92 

17-53 





11.61 

11.89 

12.25 

12.69 





75 

16.32 

16.55 

16.92 

17-44 





11.91 

12-12 

12.40 

12.76 




55-53 

80 

16. 42 

16.59 

16.89 

17.34 




66-73 

12.17 

12.31 

12-52 

12.80 




44.71 

65 

16.48 

16. Go 

16.64 

17.21 




56.11 

12.35 

12.43 

12.57 

12.76 




37-59 

90 

16.51 

16.56 

16.74 

17.05 



71-84 

48.60 

12.41 

12.44 

12.52 

12.67 



48.02 

32-55 

95 

l6.48 

16.49 

16.61 

16.85 



59-59 

43.00 

12.35 

12-32 

12-35 

12.44 



39-82 

28.78 

100 

16.42 

16-37 

16.43 

16.62 


77-17 

51.08 

38.65 

12.17 

12.09 

12.06 

12.10 


51-51 

34.12 

25.85 

105 

16.32 

16.21 

16.22 

16-35 


63.09 

44.81 

35.16 

11.91 

11-77 

11.68 

II.67 


42.09 

29-90 

23-48 

110 

16.19 

16.02 

15-98 

16.06 


53-49 

39-98 

32.29 

11.61 

11-40 

11.26 

11.19 


35.67 

26.65 

21.51 

115 

16. 04 

15.82 

15-72 

15-76 

66.60 

46-51 

36-14 

29-87 

11.31 

11.04 

10.84 

10-72 

44.41 

30.98 

24.04 

19-84 

120 

15-89 

15.61 

15-46 

15-44 

55-81 

41.19 

32-99 

27.79 

11.04 

10.71 

10.47 

10-31 

37.19 

27**40 

21.90 

18.39 

125 

15.74 

15-40 

15-20 

15-14 

48.09 

36.99 

30.36 

25.98 

10.80 

10.43 

10.15 

9.96 

32.01 

24-57 

20.10 

17.13 

130 

15-60 

15-20 

44.95 

14.84 

42.27 

33-58 

28.11 

24.38 

10.60 

10-19 

9.88 

9.67 

28.10 

22.26 

18.56 

16-02 

135 

15-46 

15.02 

14.72 

14.56 

37-71 

30.74 

26.16 

22.95 

10.44 

9-99 

9.66 

9-43 

25.03 

20-33 

17.22 

15.03 

140 

15-35 

14.84 

14.50 

14.30 

34.04 

28.33 

24.45 

21.67 

IO.31 

9.83 

9.47 

9-22 

22.56 

18.70 

16.06 

14.16 

145 

15.24 

14.69 

14-30 

14.05 

31.00 

26.25 

22.93 

20-51 

10.20 

9.69 

9.31 

9.05 

20.51 

17.29 

15.03 

13.37 

150 

15.16 

14.54 

14.11 

13. 82 

28.44 

24.44 

21.58 

19.45 

10.11 

9-58 

9-17 

8.89 

18.78 

16.07 

14.13 

12.67 

155 

15.08 

14.41 

13-93 

13.61 

26.24 

22.84 

20-35 

18.48 

10.04 

9-48 

9-05 

8-75 

17.31 

15-00 

13.31 

12.04 

160 

15.03 

14.29 

13-77 

13-41 

24.33 

21.41 

19.25 

17-59 

9-99 

9-39 

8.94 

8.62 

16.03 

14.05 

12.58 

11.46 

165 

15-59 . 

14.46 

13-61 

13.21 

22.65 

20-13 

10.24 

16.77 

10.26 

9.49 

8.88 

8.50 

H.91 

13.21 

11.92 

10-93 

170 

16.51 

15.23 

14.22 

13-42 

21.16 

18.98 

17.31 

16.01 

10.87 

10.00 

9-31 

8.76 

13-93 

12.45 

11.32 

10.44 

175 

17.52 

16.06 

14.92 

14.01 

19.82 

17.92 

16.45 

15.30 

11-53 

10.54 

9.76 

9-15 

13-04 

11-75 

10.76 

9.98 

180 

18.62 

16-95 

15-66 

14.64 

18.62 

16.95 

15.66 

14.64 

12.25 

11.12 

10.24 

9-55 

12.25 

11.12 

10.24 

9-55 


l8 




TABLE I.- COORDINATES OF THE TRACE OF THE BOUNDARY OF THE GEOMAGNETIC FIELD IN 
THE MERIDIAN PLANE CONTAINING THE DIPOLE AXIS AND THE SUN-EARTH LINE FOR A 
RING CURRENT LOCATED AT A GEOCENTRIC DISTANCE OF 60,000 KM - Concluded 
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TABLE II.- COORDINATES OF THE TRACE OF THE BOUNDARY OF THE GEOMAGNETIC FIELD IN 
THE MERIDIAN PLANE CONTAINING THE DIPOLE AXIS AND THE SUN-EARTH LINE FOR A 
RING CURRENT LOCATED AT A GEOCENTRIC DISTANCE OF 50,000 KM 


(a) Stream = 10 16 protons/cm sec 2 ; ring current = 2 million amperes 


(b) Stream = 2.5xl0 ls protons/ 


11*5 230 34.5 


11.5 23.0 34.5 



(c) Stream = 10 15 protons/cm sec 2 ; ring current = 5 million amperes 


18.10 

20.16 

22.94 

19.27 

21.66 

24.95 

20.56 

23.35 

27.28 

21.99 

25-28 

30.02 

23.61 

27-50 

33.30 

25-44 

30.11 

37-32 

27-55 

33-21 

42.37 

30.01 

36.99 

48.97 

32-92 

41.71 

57-95 

36.45 

47.80 

70.99 

1*0.83 

55-99 


46.42 

53-66 

64.25 

67.64 

87-63 

69-21 

57-42 


74.41 

49.22 


60.84 

43-19 


51.60 

38.56 

64.25 

44.68 

34.87 

53-86 

39-77 

31.86 

46.42 

35.74 

29.34 

40.83 

32.47 

27-20 

36.45 

29.75 

25-34 

32.92 

27-44 

23.72 

30.01 

25-46 

22.28 

27.55 

23.72 

20.99 

25.44 

22.19 

19.82 

23.61 

20.83 

18.76 

21.99 

19.60 

17.79 

20.56 

18.48 

16.90 

19-27 

17-46 

16.07 

18.10 

16.52 

15-29 













































TABLE II.- COORDINATES OF THE TRACE OF THE BOUNDARY OF THE GEOMAGNETIC FIELD IN 
THE MERIDIAN PLANE CONTAINING THE DIPOLE AXIS AND THE SUN-EARTH LINE FOR A 
RING CURRENT LOCATED AT A GEOCENTRIC DISTANCE OF 50,000 KM - Concluded 


(e) Stream - 2- 5xl0 16 protons/cn sec 2 ; ring current - 5 million amperes 


(f ) Stream = 10 17 protons/cn sec 2 ; ring current - 5 million amperes 



13.06 

15.47 

10.23 I 

H.43 

13-06 

15.47 

12.06 1 

14.09 


12.29 1 

14.21 

17.10 

11.17 1 

12.92 

nRJ 

13-26 1 

15.55 

19.08 

10.37 1 

11.92 


14.36 

17-12 

21.53 

9.75 ] 

11.12 

13-37 1 

15.63 

19-01 

24.64 

9.68 1 

IO.69 

14.42 1 

17-13 1 

21-32 

28.74 

9-62 

10.55 

15.64 1 

18.92 . 

24.24 1 

34-42 

9-59 

IO.43 

17.06 

21.10 

28.03 

42.83 

9-57 1 

10.34 

18.74 1 

23.82 1 

33-21 1 


9-58 

10-27 

20.79 

27-33 

40.70 


9-60 

IO.23 

23-31 

32-05 



9-65 

10.22 | 

26.54 

38.74 



9-72 

10.22 1 

30.82 




9-82 

10.25 ' 

36. eo 




9-93 

10.29 | 

45.75 




10.05 

10-34 . 





10.14 

10-37 1 




36.93 

10-19 

10-35 j 




31-05 

10.15 

10.26 . 



39-72 

26.89 

10.01 

10.07 1 



32.95 

23.78 

9-76 

9-78 i 


42.65 

28.23 

21.36 

9.44 

9-42 


34.86 

24.76 

19-42 

9-10 

9-03 

45.75 

29-55 

22.07 

17.80 

8.78 

8.69 

36.80 

25.68 

19-93 

16.44 

8-52 

8.39 

30.82 

22.73 

18.17 

15-27 

8.30 

8.15 

26.54 

20-39 

16.69 

14.24 

8.12 

7-96 

23-31 

18.49 

15.44 

13-34 

7.97 

7-6o 

20.79 

I6.91 

14-35 

12.54 

7-84 

7-66 

18-74 

15-57 

13-40 

H.83 

7.73 

7-53 

17.06 

14-41 

12.56 

11.20 

7.63 

7-42 

15-64 

13.41 

11.82 

10.62 

7-55 

7-32 

14.42 

12-53 

11.15 

10.10 

7.46 

I 7-22 

13-37 

11. 75 

10.55 

9-63 

7-44 

7-13 

12.45 

12.05 

10.00 

9-19 

7.81 

7-37 

H.63 

10.42 

9-50 

8.78 

8.19 

1 7-69 

10.89 

9-84 

9.03 

8.40 

8.60 

8.04 

IO.23 

9-31 

8.60 

8-04 

1 


(g) Stream - 10 16 protons/cn sec 2 ; ring current = 10 million anperes 


0 

12-99 

14-57 


5 

12-23 

13.60 

Bjvfijjj] 

10 

n.53 

12-73 

■hBi 

15 

10.89 

n.94 

13.29 

20 

10.59 

II.36 

12-40 

25 

10.65 

n.36 

12-30 

30 

10.72 

n.38 

12.24 

35 

10.82 

11.42 

12-21 

40 

10.93 

n.48 

12.19 

45 

n.06 

11.55 

12-20 

50 

11.20 

n .65 

12.23 

55 

n.37 

H.76 

12-28 

60 

11.54 

n.ee 

12.34 

65 

n.73 

12.01 

12.40 

70 

11.91 

12.14 

12.47 

75 

12.07 

12.25 

12-52 

eo 

12.20 

12-33 

12.54 

85 

12.28 

12-37 

12.53 

90 

12.31 

12.35 

12.46 

95 

12.28 

12-28 

12-34 

100 

12.20 

12.15 

12-17 

105 

12.07 

n.97 

U.95 

no 

11.91 

n.76 

11.70 

115 

11-73 

n-53 

11-42 

120 

n.54 

n.30 

U.15 

125 

11-37 

n.08 

IO.89 

130 

11.20 

10.87 

10.64 

135 

11.06 

10.68 

10.41 

140 

10-93 

10.51 

10-21 

145 

10.82 

10.36 

IO.O3 

150 

10.72 

10.23 

9-86 

155 

10.65 

10. n 

9-72 

160 

10-59 

10.01 

9-58 

165 

10.89 

10.07 

9-46 

170 

n.53 

10.60 

9-86 

175 

12.23 

n.17 

IO.34 

160 

12.99 

n.79 

10.85 




23.O 

34-5 

s 2 

10.19 

12.23 

11.09 

13-51 

12.13 

15-06 

13-36 

16.98 

14.84 

19-43 

16.65 

22.67 

18.94 

27-14 

21.92 

33*78 

25.98 

44.67 

31.86 




(h) Stream = 2-5xlO 10 protons/cn sec 2 ; ring current = 10 million amperes 


10.95 


14.13 

16.86 

10.30 

pifi 

13-04 

15.38 

9-72 


12.09 

14.12 

9-18 

10.08 

11-25 

13-06 

8-97 

9.69 

IO.63 

12.21 

9-03 

9*70 

10-56 

11.71 

9-10 

9-72 

10.52 

11.56 

9.20 

9-77 

10-50 

12.45 

9-31 

9-84 

10-51 

11-37 

9-44 

9-93 

10-54 

n.32 

9-60 

10.04 

10.60 

11.29 

9-79 

10-18 

10-68 

n .29 

10.00 

IO.34 

10-78 

H.32 

10-23 

10-52 

IO.89 

H- 35 

10.47 

10-70 

11.00 

11.40 

10.69 

10.87 

11.11 

11.43 

10.88 

12-00 

n.ie 

n .44 

12.00 

11.07 

11-20 

12-40 

12.04 

11.07 

u.15 

II.30 

21-00 

10.98 

n.02 

12.12 

10.88 

10.82 

10.81 

10.87 

10.69 

10.58 

10.53 

10.55 

10.47 

IO.31 

10.21 

10.19 

10.23 

10-02 

9.88 

9-81 

10.00 

9-74 

9-56 

9-45 

9-79 

9-48 

9-26 

9-12 

9-60 

9-26 

9.00 

8.83 

9-44 

9-06 

8.78 

8.58 

9-31 

8.89 

8.58 

8-37 

9-20 

6-75 

8-41 

8.18 

9-10 

8.63 

8.27 

8.02 

9-03 

8-52 

8.14 

7-87 

8-97 

8.43 

8.03 

7-74 

9-18 

8-47 

7-93 

7-62 

9-72 

8.92 

8.28 

1 7-78 

10.30 

9-40 

8.68 

8.12 

10-95 

9- 92 

9*11 

1 8.48 
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TABtE III.- 


COORDINATES OF THE TRACE OF THE BOUNDARY OF 
THE GEOMAGNETIC EQUATORIAL PLANE; A 


THE GEOMAGNETIC FIELD 

= 0 ° 


IN 



(a) Geocentric distance to ring current = 

60,000 Km 



(b) Geocentric distance to ring current » 

50,000 Km | 

\ n v 2 

\ 
deg \ 

IO 10 

2.5*10 l ° 

10 13 

10 l ° 

2.5xl0 l ° 

10 17 

10 10 

2.5xl0 LO 

10 LO 

— 
2. 5x10 10 

10 13 

10 10 

2. 5x10 10 

10 17 

IO 10 

2.5XIO 10 I 

2 million amps 

5 million amps 

10 million amps 1 

2 million ompc 

5 million amps 

10 million amps | 

90.0 

11.22 

10.49 

16.51 

12.41 

11.42 

10.46 

13-76 

12.51 

10.15 

9.20 

15.36 

11.16 

10.07 

8.90 

12.31 

11.04 

95*0 

11.22 

10.49 

16.52 

12.42 

11.43 

10-46 

13.77 

12.51 

10.15 

9. a 

15.37 

11.16 

10.07 

0.98 

12.32 

11.05 

100.0 

11.24 

10.50 

16.55 

12.44 

11.44 

10.47 

I3.8O 

12.53 

10.17 

9.22 

15.41 

11.19 

10.09 

8. 99 

12.34 

11.07 

105.0 

11.27 

10.52 

16.61 

12.48 

U.47 

10.49 

13-04 

12-57 

10.21 

9.25 

15-46 

11.22 

10.12 

9-01 

12.39 

11.10 

110.0 

U.31 

10.55 

16.70 

12.53 

11.51 

10.51 

13-90 

12.62 

10.25 

9.28 

15-55 

11.27 

10.16 

9-04 

12.44 

11.15 

115.0 

11.36 

10.59 

16.81 

12.59 

U.56 

10.54 

13-98 

12.68 

IO.3I 

9-33 

15.65 

11.34 

10.21 

9.07 

12.52 

11.21 

120.0 

11.43 

10.64 

16.94 

12.60 

11.62 

10.58 

14.08 

12.76 

IO.39 

9-3 8 

15.79 

11.42 

10.28 

9 U 

12.61 

11.29 

125.0 

U-52 

10.70 

17.11 

12.78 

11.70 

IO.63 

14.20 

12.86 

10.48 

9.45 

15.95 

11-53 

IO.36 

9-17 

12.73 

U.38 

130.0 

11.62 

10.77 

17.30 

12.90 

11.80 

10.69 

14.34 

12.97 

10.59 

9-54 

16.13 

11.65 

10.46 

9.23 

12.87 

11.50 

135.0 

11.74 

10.85 

17-53 

13.04 

11.91 

10.77 

14.50 

13-U 

10.72 

9.64 

16.36 

11-79 

10-57 

9-31 

13-02 

11.6k 

140.0 

11.88 

10.96 

17.79 

13.21 

12.04 

10.85 

14.70 

13-2? 

IO.87 

9-76 

16.61 

11.96 

10.71 

9-40 

13.21 

11.78 

145.0 

12.04 

11.08 

18.10 

13.40 

12.19 

10-95 

14.92 

13.45 

11.05 

9.89 

16.90 

12.15 

10.86 

9.50 

13-42 

11.96 

150.0 

12.23 

11.22 

18.44 

13.62 

12-37 

11.08 

15.17 

13-66 

11.25 

10.05 

17-24 

12.36 

11.04 

9.63 

13-67 

12.16 

155*0 

12.45 

11.39 

18.84 

13-87 

12.58 

11.22 

15-47 

13-91 

11.48 

10.24 

17.62 

12.61 

U.25 

9.78 

13.94 

12.40 

160.0 

12.71 

n.59 

19.28 

14.16 

12.82 

11.39 

15.80 

14.19 

H.74 

10.45 

18.05 

12.90 

11.49 

9-95 

14.26 

12.66 

165.0 

13-00 

11.83 

19.60 

14.50 

13.09 

11.59 

16.19 

14.51 

12.04 

10.70 

18.54 

13.23 

11.76 

10.15 

14.62 

12. $7 

170.0 

13-34 

12.10 

20. 38 

14.88 

13.41 

11.83 

16.63 

14.88 

12.39 

10.98 

19.10 

13.60 

12.07 

10,38 

15.04 

13.32 

175.0 

13.73 

12.42 

a. o4 

15.32 

13.78 

12.11 

17.13 

15.31 

12.78 

U.31 

19-73 

14.03 

12.43 

10.66 

15.51 

13.72 

180.0 

14.19 

12.80 

a. 80 

15.83 

14. a 

12.44 

17.71 

15.80 

13.23 

11.69 

20.46 

14.51 

12.85 

10.97 

16.05 

14.19 

185.0 

14.71 

13.24 

22.67 

16.42 

14.71 

12.83 

18-37 

16-37 

13-75 

12.13 

21.29 

15-08 

13-33 

11.35 

16.67 

14.72 

190.0 

15.32 

13.76 

23.67 

17-10 

15.30 

13.29 

19-14 

17.04 

14.36 

12.64 

22.24 

15.73 

13.88 

U.79 

17.39 

15.34 

195-0 

16.03 

14.37 

24.83 

17.90 

15.98 

13*84 

20.04 

17.81 

15.05 

13-24 

23-35 

16.48 

14.53 

12.30 

18.23 

16.06 

200.0 

16.87 

15.09 

26.18 

18.83 

16.79 

14.49 

21.09 

18.72 

15.87 

13-94 

24.63 

17.37 

15.29 

12.91 

19.20 

16.90 

205.0 

17.86 

15.94 

27.77 

19.93 

17.74 

15.27 

22.33 

19.80 

16.82 

14.76 

26.14 

18. 41 

16.19 

13.63 

20.35 

17.89 

ao.o 

19.04 

16.97 

29.65 

21.24 

18.89 

16.21 

23.80 

21.08 

17-96 

15-75 

27.93 

19-64 

17.26 

14.50 

21.71 

19- 07 

a5.o 

20.46 

18. a 

31.91 

22.82 

20.27 

17.35 

25.58 

22.63 

19-33 

16 ’93 

30.08 

a .13 

18.55 

15.55 

2305 

20.50 

220.0 

22.20 

19.73 

34.67 

24.75 

a. 96 

18.75 

27.75 

24.53 

21.00 

18.38 

32-69 

22.94 

20.13 

16.84 

25.35 

22.23 

225.0 

24.36 

a.62 

38.08 

27-15 

24.06 

20.51 

30.44 

26.88 

23. C 6 

20-17 

35-93 

25.18 

22.08 

18.45 

27.83 

24-39 

230.0 

27.09 

24.02 

42.39 

30.18 

26.73 

22.75 

33-04 

29.87 

25.67 

22.44 

40.01 

20.02 

24.55 

20.49 

30.96 

27.12 

235.0 

30.64 

27.15 

48.00 

34.13 

30.21 

25.67 

38.28 

33-76 

29.06 

25-39 

45-32 

31-71 

27-77 

23.15 

35*04 

30.68 

240.0 

35.43 

31.37 

55.54 

39-46 

34.90 

2 9-62 

44.25 

39.01 

33.63 

29.36 

52.45 

36.68 

32.11 

26.74 

40.52 

35.47 

245.0 

42.19 

37.34 

66.18 

46.98 

41.53 

35.21 

52.69 

46.43 

40.07 

34.98 

62.51 

43.70 

38.24 

31.81 

48.2? 

42.23 

250.0 

52.41 

46.36 

82.24 

58.35 

51.56 

43.68 

65.44 

57.64 

49-79 

43.46 

77.70 

54.29 

47.50 

39.49 

59.97 
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Ring current at 60,000 km 


Ring current at 50,000 km 
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Figure 1.- View of coordinate system. 
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(a) Upper signs. 

(b) Lower signs. 


Figure 2 .- Integral curves for equation ( 15 ) for a = 1, p = 1, A = 0. 




Figure 3.- Trace of the boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line: a = 1 . u, = 1 . A = 0. 
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Figure 5*- Trace of the boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line*, a = 1, u = 1, A = 11.5°. 
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Figure 6.- Integral curves for equation ( 15 ) for a = 1, jj. = 1, A = 34-5°. 
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Figure 7*- Trace of the boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line; a = 1 ; \i = 1, ~K = 3^ • 5 • 






Figure 9*- Traces of the boundary of the geomagnetic field in the equatorial plane and in the meridian 
plane containing the dipole axis and the sun-earth line, a = 1, |j. = 1, A = 0. 
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Figure 10.- Traces of the boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line for various strengths of the 
ring current; nv 2 = 10 16 protons/cm sec 2 , a = 6xl0 9 cm = 60,000 km. 
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Figure 11.- Traces of the boundary of the geomagnetic field in the equatorial plane and in the meridian 
plane containing the dipole axis and the sun-earth line for various strengths of the ring current; 
nv 2 = 10 1S protons/cm sec 2 , a = 6xl0 9 cm = 60,000 km, ?\ = 0. 




Figure 12.- Traces of the boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line for various intensities of 
the solar wind; a = 6xl0 9 cm = 60,000 km, i = 5xl0 5 e.m.u. = 5 million amperes 










Figure 13*- Traces of the boundary of the geomagnetic field in the equatorial plane and in the meridian 
plane containing the dipole axis and the sun-earth line for various strengths of the ring current; 
nv 2 = 10 16 protons/cm sec 2 , a = 6xl0 9 cm = 60.000 km, 7\ = 0. 
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Figure ]A.- Traces of the "boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line for various strengths of the 
ring current; nv 2 = 10 16 protons/cm sec 2 , a = 5xl0 9 cm = 50,000 km. 


38 






Figure 15. - Traces of the boundary of the geomagnetic field in the equatorial plane and in the meridian 
plane containing the dipole axis and the sun-earth line for various strengths of the ring current; 
nv 2 = 10 16 protons/cm sec 2 , a = 5xl0 9 cm = 50,000 km, A = 0. 
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Figure l6.- Traces of the boundary of the geomagnetic field in the meridian plane 
containing the dipole axis and the sun-earth line for various intensities of 
the solar wind; a = 5xl0 9 cm = 50,000 km, i = 5xl0 5 e.m.u. = 5 million 
amperes . 
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Figure 17 • - Traces of the boundary of the geomagnetic field in the equatorial plane and in the meridian 
plane containing the dipole axis and the sun-earth line for various intensities of the solar wind; 
a = 5xl0 9 cm = 50,000 km, i = 5xl0 5 e.m.u. = 5 million amperes, A = 0. 
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